Topologically decoherence-protected qubits with trapped ions 
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We show that trapped ions can be used to simulate a highly symmetrical Hamiltonian with 
eingenstates naturally protected against local sources of decoherence. This Hamiltonian involves 
long range coupling between particles and provides a more efficient protection than nearest neighbor 
models discussed in previous works. Our results open the perspective of experimentally realizing in 
controlled atomic systems, complex entangled states with decoherence times up to nine orders of 
magnitude longer than isolated quantum systems. 



It is universally accepted that quantum computation 
would be able to solve diverse classes of hard problems 
more efficiently than its classical counterpart How- 
ever, its practical realization is made difficult by the 
conflicting requirements imposed by the absence of de- 
coherence, qubit manipulation and scaling [3]. Gener- 
ally, larger dimensional systems are in general more sen- 
sitive to decoherence, so that scaling of quantum infor- 
mation processors is a huge experimental challenge. In 
the present work, we show how to avoid decoherence and 
protect complex multiparticle quantum states (that can 
define a qubit) from its effects in an experimentally real- 
izable trapped ion system. 

There are two conceptually different approaches to 
suppress the decoherence to a level required for quantum 
error correction. One is to reduce directly the physical 
noise that leads to decoherence. It is however very diffi- 
cult due to the stringent requirements imposed by quan- 
tum error correction [4(. Another is to encode qubits in 
decoherence free subspaces (DFS) which are decoupled 
(at least to first order) to the dominant sources of noise. 
The efficiency of the latter strategy is demonstrated by 
experiments showing a dramatic increase in the decoher- 
ence times in diverse physical systems: atom clouds [5j, 
trapped ions [6] or superconducting circuits [7] in which 
noise mostly originates in the fluctuating electric or mag- 
netic fields. In trapped ion systems, DFS are usually hy- 
perfine states that can either be symmetrically coupled 
to the magnetic field [8], suppressing dephasing to first 
order, or have the same energy [9], avoiding spontaneous 
emission. 

Here we propose a design that provides higher de gree 
decoupling from environment. The main idea [ill. Il2l. 
3 EI EE ES E3 of our approach is to implement a 

long range Hamiltonian which non-local symmetries en- 
sure that not only the linear but also quadratic, cubic, 
etc coupling to the environment of the lowest doublet 
vanishes, providing thereby a high degree of protection 
against the effects of most sources of physical noise. The 



Hamiltonian introduced here is more efficient for deco- 
herence protection than previously studied ones, as is 
shown below. In our scheme, the protected states are 
complex entangled states of N 2 spin 1/2 like particles. It 
is an example of multidimensional quantum state that is 
far more robust against decoherence than its individual 
elements. 

In the second part of this letter we propose the im- 
plementation of the model Hamiltonian in trapped ion 
systems. The long range interaction responsible for pro- 
tection in our Hamiltonian is particularly adapted to 
trapped ions systems. We show how this type of experi- 
mental system that has been already employed to imple- 
ment basic operations of quantum logics [10[ can be also 
used to realize the protection by topological order. 

We start by introducing the long range Hamiltonian 
acting in a collection of two level systems. 
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Here we label each particle by its position in a two di- 
mensional lattice: a* J are Pauli matrices of the spin sit- 
uated at the intersection of the i th row and j th column. 
This Hamiltonian couples the particles in the same row 
by afjaf k interaction and those in the same column by 
a i ' k a j k interaction. Consider the projections of all spins 
in a given row onto ?/-axis. Evidently the projection is 
preserved by the column interaction while row interaction 
might change it by or 2 units. Thus all terms in the 
Hamiltonian preserve the parity of the projection of 
the spins in one row. Mathematically, it can be described 
by non-local symmetry transformations generated by the 
operators Pi = a ij an d Qj = Yli crfj that involve 
product of all spins in each row or column. It can be eas- 
ily shown that the operators of a given set commute with 
all others in the same set, i.e. [Pi, Pj] = = [Qi, Qj] but 
anti-commute with all the operators of another sets, i.e. 
PiQj + QjPi = 0. They also satisfy if = 1 = Q), Vi, j 
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and [H, Pi(Qj)] = 0. These commutation relations imply 
the existence of doubly degenerate states, combining the 
advantages of previously introduced DFS 0, @, @]. In- 
deed, a state corresponding to an even value of the pro- 
jection along the ?/-axis along one row is converted into 
an odd parity state by the action of Qj operator that 
commutes with the Hamiltonian. Finally, the numer- 
ical diagonalization of the Hamiltonian (pQ) shows that 
the two ground states are separated from the rest of the 
spectrum by a gap A which depends weakly on the size 
of the system. 

We now discuss the effect of a generic local noise de- 
scribed by the operator 
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where the b i j are arbitrary time-dependent coefficients. 
Note that no other hypothesis is needed concerning the 
noise besides its local character: this noise can be due 
to energy fluctuations caused by random fields acting on 
the atomic system (described by a z operators) as well 
as spontaneous emission (described by <r_ operators). 
Equation ([2]) should be added to the model Hamiltonian 
(PQ). Individual terms of ([2]) do not commute with some 
of the symmetries, Pi and Qj, partially lifting the ground 
state degeneracy. When only k < N such terms are 
present 2N — 2k symmetries remain intact ensuring the 
degeneracy of the ground states levels. Thus, decoher- 
ence is due to the simultaneous action of N noise terms, 
in a processes taking place in the N th order perturbation 
theory. The small parameter controlling the perturba- 
tion development is max(6^' z )/A. Decoherence rate ac- 
quires an additional small factor of (max(6^ 2y ' 2 )/A) iV 
when N noise terms act simultaneously. Note that the 
decoherence rate of a system consisting of N 2 interacting 
particles usually increases with the number of particles, 
a behavior which is in clear contrast to the one observed 
here. The resulting decoherence rate is thus 



a N T (max(&?f*)/A; 



N-l 
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where is a numerical constant of the order of the 
unity. For this protection to be efficient the gap A should 
remain large as N increases. 

The results of the exact numerical diagonalization of 
the long range Hamiltonian (pQ) with J x — J y reported in 
Tab. [J show that this is indeed the case. This property 
of (PQ) is a clear advantage over short range models con- 
sidered in lH, 3, 15, 16j]. It implies that the suppression 
of the noise is more efficient in the present system. In- 
deed, numerical calculations show that the static fields 
b*j ,z randomly distributed in the interval (—0.1, 0.1) J 
result in a very small splitting of the degenerate doublet 
SE ~ 10 -4 — 10 _5 J , two orders of magnitude better 
than that of [14[ in agreement with the estimates above. 



We now discuss possible implementations of the math- 
ematical model JT) in a physical system. As was shown 
previously, short range Hamiltonians with topologically 
protected doublets can be realized by Josephson junction 
arrays [H,[l3- Although promising, this approach is dif- 
ficult due to the required high degree of similarity of a 
large number of nanoscale junctions. Very recently it was 
shown that it is possible to realize these Hamiltonians by 
polar molecules in optical lattices, [l6[ but this is also 
technologically difficult. The long range Hamiltonian (pQ) 
has the advantage, besides the stronger protection de- 
scribed above, that it can be realized in a trapped ion 
system with present day technology using laser induced 
coupling proposed in [18[ (see Fig. [1]) . This type of in- 
teraction has already been used to create multiparticle 
entangled states in [211 ] . 

In this approach the physical system consists of N 2 
ions that can either be arranged in a two dimensional 
array of microtraps, see Fig. [TJ or in a linear trap, see 
Fig. [2j In both situations, strings of ions are illuminated 
by laser fields tuned close to a particular resonance fre- 
quency of individual ions|18L Il9[ l20j|. In these conditions 
the relevant degrees of freedom of each ion can be repre- 
sented by a pseudospin 5=1/2 variable. As discussed 
below one can choose the laser fields so to generate an 
effective Hamiltonian of these pseudospins coupling ions 
two-by-two. Such a pseudospin can be implemented us- 
ing an electronic ground state and an excited level di- 
rectly coupled by a single photon transition (22[, in this 
case spontaneous emission will not bring the system out 
of the pseudo spin space. Furthermore, for long lived 
states of alkaline earth ions [28], the corresponding weak 
metastability induces a small noise term of the form ([2]), 

///// 




FIG. 1: Two laser beams of slightly different frequencies (blue 
and red) produce virtual processes that change the state of 
one ion and create (annihilate) a phonon mode of N x N lat- 
tice. Each of these laser fields alone does not produce a real 
ion transition (with or without a phonon mode emission or 
absorption). The sum of their frequencies is exactly equal 
to the energy needed to change the state of two ions simul- 
taneously which is described by the cr^^a^^ in the effective 
Hamiltonian. If the phonon mode involved in this process 
corresponds to global translational motion of all ions (center- 
of-mass mode), the effective pair- wise interaction does not 
depend on the distance between two ions. It can be restricted 
to the desired row (column) form if only a subset of ions is 
illuminated at a given time. The effect of the sequential ap- 
plication of the beams to rows and columns is described by 
the effective Hamiltonian [T] 
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which is efficiently treated by the incipient protection 
mechanism of Hamiltonian (pQ). 

We now sketch the derivation of the effective Hamil- 
tonian (p] ) and estimate realistic constraints on its pa- 
rameters. Each ion is in one of two quantum states 
with an energy difference uo egi the interaction between 
ions produce a collective vibrational mode of frequency 
v corresponding to the global displacement of all ions. 
The ions are subjected to two laser fields of frequencies 

^1,2 



± (y + 8) and Rabi frequency Q. In the inter- 
acting representation the physical Hamiltonian describ- 



ing the ion system in the to 
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+ [y + 5) laser field is 
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h.c. 

(4) 

where rj is the Lamb-Dicke parameter defined as 
(h 2 k]/2Mhv) 1 / 2 and a and are the annihilation and 
creation operators of the phonon mode. Provided that 
Q <C v and r) <C 1 the main contribution to the effec- 
tive Hamiltonian comes from the terms oscillating with 
frequency S. For small interaction constant 77^ <C S 
these terms can be treated in the perturbation theory 
resulting in the desired x( J ' n ) 2 term of the effective 
interaction where n is a unit vector whose direction in 
xy plane is controlled by the relative phases of the laser 
fields [18]. Generally, a stronger interaction might lead to 
the excitation of the phonon mode. This process can be 
suppressed [19] by choosing the time of the interaction r 
in such a way that the phonon system returns to its ini- 
tial state, t5 = 2ttK. In this case the effective interaction 
remains x(J • n) 2 with \ — 

Because the interaction occurs only between the ions 
illuminated at one time, the spatial form of the physical 
array does not need to be directly related to the col- 
umn/row form of the effective Hamiltonian (pp). Thus, as 
we illustrate in Fig. [21 a possible configuration rendering 
the implementation of the Hamiltonian easier in existing 
systems, consists of using a linear array of ions. 



2x2 



3x3 



4x4 



5x5 



SRI 
LRI 



0.84 J x 
0.84 J x 



0.58 J x 
0.96 J x 



0.32 J x 
0.92 J x 



0.20 J x 
0.80 J x 



TABLE I: Gap as a function of the array size for short range 
interaction (SRI) and long range interaction (LRI) 

This model also assumes that the lasers interact only 
with mode of frequency v. This becomes difficult to 
achieve in large systems where the number of modes is 
large. Small coupling to these modes will induce a po- 
sition dependent coupling in the effective Hamiltonian 
modifying only slightly the gap in the spectrum. Large 
coupling to these modes would lead to their excitation 
which is much more dangerous. To avoid this process we 
need to ensure that max(J, Q) < v — v\ where v\ is the 
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FIG. 2: An array of 3 x 3 implemented by a linear trap with 
9 ions. To generate the interaction in row 1 (insert a) one 
applies laser light to ions 1, 2, 3 (insert b) while column 2 
interaction requires light on ions 2, 5 and 8. (insert c). The 
advantage of this scheme is that it does not require micro 
trap fabrication. However, it becomes difficult to implement 
for large arrays (N > 3) because the distance between the ions 
in the center of the array gets smaller making their individual 
addressing difficult. 



frequency of the closest mode. In the case of a string 
of ions, v — v\ « (1 — y/3)v and is roughly independent 
of the number of ions [23], for a fixed trapping potential. 
Thus, in order to maximize interaction, we need to adjust 
the parameters so that v — v\ is maximal while keeping 
the ions at distances sufficiently far apart so that they 
can be addressed individually. This leads to conditions 
on the translational mode frequency, which can be satis- 
fied experimentally even for chains of 8 ions j24[. Using 
the results of 23|, [25[ for the phonon spectrum of a sin- 
gle ion chain and assuming a fixed inter- ion distance, we 
find that v — ^l^^-xlO MHz, which implies that such 
implementation is feasible for N < 3. 

To implement larger arrays one can use a two dimen- 
sional configuration of surface traps of a few meV depth, 
similar to those used in [26]. When the traps are at 
a distance a m i n from each other the Coulomb repulsion 
becomes of the same order of magnitude as the trapping 
potential. In this case, one can achieve a frequency of the 
translational modes of the order of v = 10 MHz [26[ and 
our numerical results show that for a 5 x 5 square array 
configurations v — v\ « 0.1 Coupling between rows and 
columns is then avoided by applying lasers that alternate 
between lines and columns. The resulting Hamiltonian 
is the sum of the terms in (pQ) provided that the dura- 
tion of each pulse satisfies rJ XyV < 1. We summarize the 
experimental values and resulting couplings in Tab. [Til 



4 ions 



9 ions 



5x5 ions 



tt(Hz) 
6 (Hz) 

J X ,y(Hz) 



io 4 

10 4 



io 4 
io 4 



io 5 
io 5 



TABLE II: Parameters and induced coupling for the various 
arrays, v — 1 MHz for 4 and 9 ions and 10 MHz pour 5x5 
ions, 77 = 0.1 in both cases, J x , y — r] 2 Q 2 /S, K = 1. 



The strength of the induced interaction, J x ,J y and 
thus the gap in the spectrum A of the Hamiltonian (pQ) 
is limited by the condition that the laser do not lead 
to the excitation of spurious phonon modes. Thus, to 
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achieve a maximal interaction strength one needs a vi- 
brational spectrum with the largest possible gap between 
the global translational mode and the rest of the spec- 
trum. Another constraint comes from the condition that 
the distance from the ions is sufficiently large to allow 
for their individual addressing by laser beams. This gap 
is very sensitive to the geometry of the lattice. We have 
identified the two most promising candidates mentioned 
above: small (4 and 9 ions) one dimensional (Fig. [2]) and 
larger (5 x 5) highly symmetric two-dimensional struc- 
tures (Fig. HI). 

We now estimate the effective decoherence rates which 
can be obtained in these three systems in realistic con- 
ditions. Assuming that the allowed minimal separation 
between ions is 2/im we get (Tab. (TTJ) that the induced 
couplings J XiV of the effective Hamiltonian are 10 4 Hz 
for the linear array and 10 5 Hz for the square array. As 
shown in Tab. H the gap is of the same order of mag- 
nitude as the coupling coefficient in the isotropic case, 
J x — Jy m 

The effective decoherence rate ((3]) is determined by 
the largest noise term and by the single ion decoherence 
rate. In an ideal system where the only origin for noise 



would be the individual ion decoherence 



To), 



the effective lifetime would be astronomical (assuming 
a minimal separation between ions a m i n = 2/xra we get 
10 25 s for a 5x 5 array). In a more realistic situation, 
one must take into account noises induced by the ions 
manipulation with laser light: the dominant noise turns 
out to come from laser frequencies noise Sf. Frequencies 
drift with time resulting in a H = HSf ■ afj term in 
the effective Hamiltonian,. Though one can suppress this 
noise down to the level of some Hz [27] we shall use a more 
realistic value of 500 Hz for our estimates summarized in 
Tab. [Till 



4 ions 



9 ions 



5x5 ions 



^eff(Hz) 

r(s) 



1.5 • 10" 

6.6 • 10 2 



7.5 • 10- 
1.3 • 10 4 



1.9 • 10 
5.3 • 10 10 



TABLE III: Effective decoherence rates and topological qubit 
lifetimes. 



Global state initialization can be achieved by first 
imposing a large effective field along the x direction, 
yielding a large additional term proportional to H ext = 
±^<rf^ in (pQ). This leads to the global state of the 
array jo) = [1 |0)»j or I 1 ) = EII 1 )^ depending on the 
overall sign in H ext . These states satisfy Qj |0) = |0) and 
Qj |1) = - |1) (if N is odd). When H ext is switched off 
adiabatically (i.e on a time scale much longer than 1/A), 
the final Hamiltonian becomes (pp). At all stages of this 
evolution the Hamiltonian commutes with Qj operators. 
Thus, the two initial unprotected states |0) and |1) evolve 
into the protected states with the same quantum num- 
bers. Note that Pi operators commute only with the final 



Hamiltonian, and the initial states are not eigenstates for 
them. The final state can be reconstructed using stan- 
dard methods [24] measuring the state of individual ions 
in the <j x or cr y basis, after the protecting interaction (pQ) 
has been switched off. 

We have shown that it is possible to use the flexibility 
of ion traps to form long living quantum states topologi- 
cally protected from decoherence, even small 4 ion strings 
being expected to give relaxation times one order of mag- 
nitude longer than currently achieved (Tab. ITTTj) . It is 
quite likely that similar techniques can be used for exper- 
imental implementation of more exotic quantum states 
in controlled atomic systems, for instance the states de- 
scribed by the anyon statistics that would make possible 
to implement a full quantum computation in the pro- 
tected subspace. 
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multiplets have a typical lifetime of 1 s 



